Creep motion of an intruder within a granular glass close to jamming 
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We investigate experimentally the dynamics of an intruder dragged by a constant force in an 
assembly of horizontally vibrated grains close to jamming. At moderate packing fractions, the 
intruder moves rapidly as soon as the force is applied. Above some threshold value of the packing 
fraction which increases with the applied force, the intruder exhibits an intermittent creep motion 
with strong fluctuations reminiscent of a "crackling noise" signal. These fluctuations behave in 
a critical manner at the jamming transition unveilled in a previous study . The transition 
separates a regime where the intruder motion is dominated by local free volume rearrangements 
from a regime where the instantaneous displacement field is strongly heterogeneous and resemble 
the force chains patterns observed in dense granular packings. 
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The understanding of mechanical properties of amor- 
phous media such as granular media, foams, emulsions, 
suspensions or structural glasses raised a formidable in- 
terest in the past decades [2. At high packing fractions, 
such materials eventually jam and sustain a finite shear 
stress before yielding [H [SJ [Q. Several experiments on 
granular systems have emphasized the role of dynamical 
heterogeneities in the glassy like increase of the struc- 
tural relaxation time [7 . More recently it was shown 
that dynamical heterogeneities, albeit of another type, 
also control the time scales of the jamming transition of 
a horizontally vibrated granular monolayer [1 . However 
these experiments don't provide any direct measurement 
of a mechanical response function. Conversely many ex- 
periments report on the stress-strain relation in dense 
granular packings [8 , but don't have an easy access to 
the dynamics of individual grains. Investigating the drag 
of an intruder, the subject of this paper, is a possible 
way of bridging the gap between dynamical and mechan- 
ical properties. Experiments in colloids [5], foams [4 
and granular media [H [TOl [TTJ [12] as well as simula- 
tions of structural glasses p!3] were performed along this 
line. For loose packings and large drag, that is in the 
so-called fluidized regime, the velocity dependence of the 
drag force F follows Stokes' law, F oc V [9 . For denser 
packings, experiments report either F ~ est or F oc ln(V) 
HllOlIIIlIig when the velocity is fixed, and F = Fy + V^, 
with 7 < 1 and Fy a finite yield force, when the force 
is imposed [H [13]. In both cases, the dynamics has 
been described as very intermittent. Stress fluctuations 
have been investigated in details in [H], but very little 
is known about the displacement field and the velocity 
fluctuations. 

In this letter we study the dynamics of an intruder 
dragged by a constant force, in a bidisperse mono-layer 
of horizontally vibrated grains. The experiment is run in 
the same set up (fig. [T]) and following the same protocol 
as in [T], where the jamming transition has been identi- 
fied without ambiguity: at the pressure measured in 
the absence of vibration vanishes and dynamical hetero- 



geneities exhibit a critical behaviour. Here, we observe 
that close to ^j, the intruder motion is strongly inter- 
mittent, with widely fluctuating velocities. We analyze 
these fluctuations which resemble very much a "crack- 
ling noise" signal [M] and show that the jamming transi- 
tion is signed by critical fluctuations of the intruder mo- 
tion. Investigating the displacements and the free volume 
fields around the intruder, we conclude that the transi- 
tion separates a regime dominated by local free volume 
rearrangements from a regime dominated by the rear- 
rangements of the force network. We also show that 
the above transition is distinct from fluidization [12 , ob- 
served at a looser packing fraction which depends on the 
applied force, when the intruder recovers a continuous 
motion. 

The experimental setup (fig.[l]) has been described else- 
where and we shall only recall here its most important 
elements together with the modifications imposed by the 
drag of the intruder. A monolayer of 8500 bi-disperse 
brass cylinders of diameters dsmaii = ^l^d^g = 4 ± 0.01mm 
lays out on a horizontal glass plate vibrated horizontally 
(/ = lOHz, A = 10mm). The grains are confined in a 
cell fixed in the laboratory frame. The packing frac- 
tion, ^, can be varied by tiny amounts {S(j)/(j) ~ 5.10"^) 
and the pressure exerted on the moving lateral wall is 
measured by a force sensor. The intruder consists in a 



larger particle of same height {d 



intruder 



2.4 



i) in- 





-5 5 10 

FIG. 1: Left: a monolayer of disks is vibrated horizontally, 
while dragging an intruder at constant force. Both the intrud- 
ers and the surrounding grains are tracked by a CCD camera; 
Right: a strongly intermittent and heterogeneous response is 
observed, see text for details 
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FIG. 2: From flow to jammed states: Left Parameter space, 
(force F - relative packing fraction e = (0- j at low cj) 
and large F,(o) the packing is fluidized, the intruder motion 
is continuous and F oc V] large cj) and low F (0), the 
intruder exhibits an intermittent motion and F oc \nV. The 
horizontal and vertical doted lines indicate the path followed 
in the parameter space in the present study. Right F versus 
V along the paths 0i(v) , 02 (o), 03 (^); Inset: same in log- 
lin. 



troduced in the system and pulled by a mass via a pul- 
ley perpendicularly to the vibration. The fishing wire, 
which stands over the other grains, doesn't disturb the 
dynamics. Most of the results presented here are re- 
lated to experiments performed at a constant force {Fi = 
0.677V, F2 = 1.487V, F3 = 2.627V) and varying the packing 
fraction, but we also conducted experiments at constant 
packing fractions, (0i = 0.8383,02 = 0.8394,03 = 0.8399) 
increasing the force from 0.677V to 14.257V. For compara- 
ison, the total weight of the grains is 23.117V and the 
force registered at the wall when the grains are highly 
compressed is of the order of 507V. The time unit is set 
to one plate oscillation while the length unit is chosen to 
be the diameter of the small particles. 

Starting from a low packing fraction 0, we gradually 
compress the system until it reaches a highly jammed 
state following the same protocol as in [1]. Then we 
stepwise decrease the volume fraction. In the absence of 
intruder, it was shown that the average relaxation time 
increases monotonically with the packing fraction, while 
the dynamics exhibits strong dynamical heterogeneities, 
the length-scale and time-scale of which exhibit a sharp 
peak at an intermediate packing fraction. The pressure 
measured at the wall in the absence of vibration falls to 
zero precisely below that packing fraction, hence called 
the jamming transition 0j. In the present study, the in- 
truder is inserted at its initial position in place of one 
big and two small grains before each downward step in 
packing fraction and the system is kept under vibration 
until the pressure has recovered its value in the absence 
of the intruder. Only then the force is applied on the 
wire and the intruder is dragged through the cell, while 
its motion together with that of a set of 1800 surrounding 
grains in the center of the sample is tracked by a digital 
video camera triggered in phase with the oscillations of 
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FIG. 3: Top- Left Intruder instantaneous displacements 
6x(t) for two different packing fractions. Middle-Left Distri- 
butions of Sx for 15 packing fractions. Inset: same in log- log. 
Bottom-Left Inverse average l//i+ (o - right axis) and stan- 
dard deviation over average <j+Iijl+ (♦ - left axis) of 5^ as a 
function of the reduced packing fraction e for the three ap- 
plied force Fi (red), F2 (blue) and F3 (green). Top- Right 
Cumulated distributions of the durations T of the displace- 
ment bursts for different packing fraction ; Inset: Cutoff of 
the distributions 0{e) versus 0. Middle- Right Cumulated 
distributions of the rescaled size L/Lo of the bursts for differ- 
ent packing fraction; Inset: Cutoff of the distributions ^(e) 
versus 0. Bottom- Right Scaling of the size L of the bursts 
with their duration T. The avalanche sizes are rescaled by a 
factor Lq. Inset: Lq and V (dark dotted line) follow the same 
evolution with 0. (In all case, F = F2 when not specified.) 



the plate. 

At low packing fraction, and large enough force, the 
intruder motion is continuous. When increasing the 
packing fraction, it becomes intermittent above some 
threshold, which increases with the applied force (fig|2]- 
a). Performing experiments at a given packing fraction 
and increasing the force, one observes that the applied 
force is proportional to the average velocity of the in- 
truder F oc y in the continuous motion regime, whereas 
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oc InV in the intermittent one (fig[2}b). We thereby 
identify this transition with the fluidization one [12^ . 
Note that for the largest packing fraction ^2 and ^3, we 
could not observe the fluidization. As a result, there is 
a very strong contrast between the continuous motion of 
the intruder observed at = ^1 for large enough forces, 
say F = ION, and the strongly intermittent one observed 
dit (j) = (j)2 for the same large force. Typically the intruder 
averaged velocity looses three orders of magnitude and its 
velocity flucutations gain more than five orders of mag- 
nitudes, while (j) is increased by less than 1%, already 
suggesting the existence of a sharp transition. 

We now focus on the experiments performed at con- 
stant and rather small force. When looking at the in- 
truder displacements 6x along the dragging direction 
during one vibration cycle (fig. [sj-top-left), one imme- 
diately notices very strong fiuctuations, with bursts 
of widely fiuctuating magnitude. More quantitatively, 
the probability density functions of (fig. [sj-middle- 
left) exhibits an important skewness towards the pos- 
itive displacements. We characterize the positive part 
of the distribution i.e. the displacements in the direc- 
tion of the drag force, 5x'^ , by computing the average 

value /i"^ = (dx^) and the relative fiuctuations cr^ / /J^^ = 
.1/2 " 



{^{6x^ - ji^YY'^ I ji^ . One observes (fig. 
that 1 //i"^ increases continuously by three 



bottom-left) 
orders of mag- 
nitude, while varying the packing fraction of only a few 
percent, Scjyjcj) = 2.10"^, and that cr'^ / /J^^ exhibits a peak 
at an intermediate packing fraction. Both behaviours are 
directly reminiscent of what has been recalled above for 
the dynamics in the absence of intruder. We could check 
indeed that the peak observed in the fiuctuations of the 
intruder motion coincides with a vanishing pressure in the 
absence of vibration and thereby locate it at (j)j without 
ambiguity. Note however that the precise value of de- 
pends on the precise packing that has been selected when 
the system has been compressed and that the small and 
non monotonous variations of 0j = 0.8369,0.8386,0.8379 
for the three forces Fi , F2 , F3 together with the small 
difference with the value (j)j = 0.842 reported in [T] must 
be attributed to differences in the initial conditions and 
more generally in the compression protocol (see [T5j [16] 
for a detailed discussion in the case of hard spheres). 
Figure [3]-bottom-right displays both l/ja'^ and cr"^//i"^ as 
a function of e = the relative distance to for 

the three dragging forces. The peak in the fiuctuations 
is clearly separated from the divergence of l/ja^ and is 
a new signature of the jamming transition independent 
from the one provided by the study of the dynamical 
heterogeneities. 

To complete the characterization of the intruder mo- 
tion, we analyse the duration T and the size L associ- 
ated with the displacement bursts, as is commonly done 
in crackling noise or Barkhausen noise experiments [M]. 
We impose a given reference level and define bursts as 
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FIG. 4: Displacements and free volume around the intruder 
(the intruders goes from left to right; F = F2). Top-Left In- 
terpolated average displacement field : the darker, the faster. 
Top-Right y-profiles of the average velocity along the drag 
direction for different packing fractions. Bottom-Left Av- 
erage free volume field (j) = Bottom-Right Average 
free volume vs e in front (red v) and behind (blue A) the 
intruder. 



zones where 5x is above this level. The duration T of a 
given burst is defined as the interval within two succes- 
sive intersections of 5^ with (5^, while the size L is defined 
as the integral of 8x between the same points. Despite 
some weakness of our statistics, the cumulated distribu- 
tions of T (fig |3]-top-right), respectively of L/Lq (fig |3]- 
middle-right), where Lo{(j)) is a scaling factor following 
the same evolution as (</>), can be described as power 
laws truncated by a scaling function : T~"/(T/6>(0)), 
resp. (L/Lo)A(^/^(^)). The cut-off l9(0), resp. ^(^), 
estimated by computing the experimental averages of 
T, resp. L/Lq display a diverging behaviour at 0j: 
0{(p) ~ \<p-(j)j\~^, resp. ^((/>) ~ Our statistics 

are not large enough to extract precisely the exponents a, 
T] and v; however the dynamical exponents defined 
by L/Lo(^) - T^l' (figlJtop-bottom) is reasonably well 
determined, z = 2/3 and estimates of a = 1/2, /3 = 1/3, 
T] = 2/3 and u = 1 are consistent with the data and satisfy 
the relations az = P and rj = vz. Also, the same anal- 
ysis performed on the kinetic energy of the surrounding 
grains - not shown here - is consistent with the above 
determination. 

Finally, we characterize the dynamics around the in- 
truder. The averaged displacement field (fig. |4]-top- left) is 
composed of two symmetric recirculation vortices. This 
pattern is very robust as evidenced by the shape in- 
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variance of the displacement profiles along the direction 
perpendicular to the drag (fig. [4]-top-right). In particu- 
lar, the exponential decay of the displacement amplitude 
keeps the same characteristic length across the transition 
(see inset). Such a smooth and continuous behaviour is 
in contrast with the existence of the sharp transition de- 
scribed above and must be related to the similar absence 
of signature of the transition when considering the av- 
erage relaxation time or Again the transition is 
to be found in the strongly heterogeneous instantaneous 
displacement field (fig. [T]-right), which exhibits charac- 
teristic chain-like motions. This tendency of forming 
chain-like motions is strongly enforced for packing frac- 
tion larger than We also compute the averaged free 
volume, extracted from the Laguerre's tessellation of the 
packing, around the intruder (figure |4]-bottom-left) and 
observe a small asymmetry between the front and the 
back of the intruder. Computing for instance the aver- 
age free volume in a small window in front of and behind 
the intruder (see figure [ij-bottom-right), one sees that be- 
low (j)j there is a significant excess of free volume leaving 
like a "wake" far behind the intruder, whereas above 
this asymmetry rapidly vanishes together with a strong 
decrease of the available free volume. 

Altogether, the following pictures emerge. The jam- 
ming transition, which is marked by a critical behaviour 
of the fiuctuations of the intruder displacements sepa- 
rate a regime where the intermittent motion is domi- 
nated by rapid reorganization of the free volume from 
one where the motion takes the form of chain-like struc- 
tures, which are very much reminiscent of the strong force 
network suggesting the dominant role of the stress fiuctu- 
ations [IT]- As compared to previous studies of the mo- 
tion of an intruder in a granular packing [9l [TOl [TTJ [12] , 
we have confirmed the transition from a linear viscous- 
like dependence of the dragging force with the velocity in 
the fiuidized regime to a logarithmic dependence in the 
intermittent one. The originality of the present study 
is to have unveil critical features inside the intermittent 
regime and associate them with the jamming transition. 
From a more general point of view, the observation of 
"crackling noise" statistics suggests a possible deeper cor- 
respondence in the underlying physical properties with 
other intermittent phenomena such as the creep motion 
observed under yield stress in amorphous media [18] , the 
subcritical material failure [19] and more generally, the 
pinning-depinning transition [20]. Further investigations 
in the case where the intruder is dragged at constant ve- 
locity through a spring, should provide new insight in this 
matter. We believe that performing the kind of micro- 
rheology experiment sketched out in the present study is 
a promising path for a better understanding of the mech- 
anisms at play at the jamming transition of frictional 
systems. 
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